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In this paper we investigate the possibility of a dynamical definition of an effective temperature 
for compacting granular media in the framework of the Fluctuation-Dissipation (FD) relations. We 
have studied two paradigmatic models for the compaction of granular media, which consider particles 
diffusing on a lattice, with either geometrical (Tetris model) or dynamical (Kob-Andersen model) 
constraints. Idealized compaction without gravity has been implemented for the Tetris model, and 
compaction with a preferential direction imposed by gravity has been studied for both models. 
In the ideal case of an homogeneous compaction, the obtained FD ratio is clearly shown to be in 
agreement with the prediction of Edwards' measure at various densities. Similar results are obtained 
with gravity only when the homogeneity of the bulk is imposed. In this case the FD ratio obtained 
dynamically for horizontal displacements and mobility and from Edwards' measure coincide. Finally, 
we propose experimental tests for the validity of the Edwards' construction through the comparison 
of various types of dynamical measurements. (PACS: 05.70.Ln, 05.20.-y, 45.70.Cc) 



I. INTRODUCTION 

Granular materials pi |^, |3| play a very important role 
in many fields of human life and industrial activities, like 
agriculture, building, chemistry, etc. Their properties are 
interesting not only for practical reasons, but also from 
the point of view of fundamental physics. In fact, in spite 
of their apparent simplicity, they display a wide variety 
of behaviours that is only partially understood in terms 
of general physical principles. 



Typical problems to face in the study of 
compact granular matter 



equations do not leave the microcanonical or any other 
known ensemble invariant. Moreover, just as in the case 
of aging glasses, the compaction dynamics j|, [|, [| does 
not approach any stationary state on experimental time 
scales (at least at small enough forcing) and these sys- 
tems exhibit aging |g, |, || and memory Q [□], fl2|. A 
granular media lives then always in non-equilibrium con- 
ditions. Even when one observes some stationary state 
as the result of a specific dynamics (energy injection) 
imposed to the system, one is never able to establish 
some sort of equipartition principle ruling how the en- 
ergy injected is redistributed among the different degrees 
of freedom of the system. 
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The common wisdom about granular materials defines 
them as non-thermal systems, since thermal energy can 
be generally neglected if compared to mechanical energy 
due to gravity and other external energy sources usually 
acting on these systems. In addition, a fundamental role 
in the dynamics is played by the mechanical energy dissi- 
pation due to friction and collisions among the grains and 
with the container walls: motion can take place only by 
continuously feeding energy into the system, that oth- 
erwise would get stuck into some metastable state, no 
longer exploring spontaneously the configuration space. 
Consequently, as a matter of fact, the dynamics of gran- 
ular matter is always a response to an external pertur- 
bation and in general the response will depend in a non- 
trivial way on the rheological properties of the medium, 
on the boundaries, on the driving procedure and, last but 
not the least, on the past history of the system. 

From the non-thermal character of these systems a lot 
of consequences can be drawn. The first one concerns the 
lack of any ergodicity principles: a granular media is not 
able to freely explore its phase space and the dynamical 



Despite all these difficulties, since granular systems in- 
volve a large number of particles, one is always natu- 
rally inclined to treat them with methods of statistical 
mechanics. In particular one of the main questions con- 
cerns the very possibility to construct a coherent ther- 
modynamics and on this point the debate is wide open. 
The construction of a thermodynamics would imply the 
identification of a suitable distribution that is left invari- 
ant by the dynamics {e.g. the microcanonical ensemble), 
and then the assumption that this distribution will be 
reached by the system, under suitable conditions of 'er- 
godicity'. As already mentioned since in granular media 
energy is lost through internal friction, and gained by 
a non-thermal source such as tapping or shearing, this 
approach seems doomed from the outset. Nevertheless 
one could ask whether some elementary thermodynamic 
quantities are well defined and what is their meaning. It 
is in this spirit that in this paper we address the question 
of the definition of an effective temperature for compact 
granular media. Before going in the details let us briefly 
review the state of the art on these subjects. 



2 



B. State of the art 

A lot of approaches fll3|-|^2t have been devised in the 
last years to provide a coherent scenario but till now the 
situation is quite uncertain mainly because at a funda- 
mental level there is no general argument showing that 
a particular construction should lead to the relevant dis- 
tribution for the dynamics (as one does in the case of 
conservative dynamics, leading to thermodynamics). 

Many models have been proposed in order to reproduce 
the rich phenomenology observed in granular compaction 
experiments [|], pj, but a general thermodynamical-like 
framework, based on the idea of describing granular ma- 
terial with a small number of parameters, is however still 
lacking. 

A very ambitious approach, aiming at such a descrip- 
tion of dense granular matter has been put forward by S. 
Edwards and co-workers JlJ, [l^] , by proposing an equiv- 
alent of the microcanonical ensemble: macroscopic quan- 
tities in a jammed situation should be obtained by a flat 
average over all blocked configurations (i.e. in which ev- 
ery grain is unable to move) of given volume, energy, etc.. 
The strong assumption is here that all blocked configura- 
tions are treated as equivalent and have the same weight 
in the measure. 

Very recently, important progresses in this direction 
have been reported in various contexts: a tool to system- 
atically construct Edwards' measure, defined as the set 
of blocked configurations of a given model, was proposed 
in [^| |24| ; it was used to show that the outcome of the 
aging dynamics of the Kob- Andersen (KA) |2^] model (a 
kinctically constrained lattice gas model) was correctly 
predicted by Edwards' measure. Moreover, the validity 
and relevance of Edwards' measure have been demon- 
strated ]26| for the Tetris model |2l[] , for one-dimensional 
phenomenological models p7f , for spin models with "tap- 
ping" dynamics |^8| , and for sheared hard spheres [^9| . 

At the time being however, the correspondence be- 
tween Edwards' distribution and long-time dynamics is 
at best checked but does not follow from any principle. It 
is therefore important to continue to explore its range of 
validity and therefore to test its applicability to various 
kinds of models. 

Another important message emerging from these stud- 
ies concerns the link between Edwards' approach and 
the outcomes of the measurements of the Fluctuation- 
Dissipation relations. In j2^, |24|, it has been shown in 
the framework of two non-mean field models, the Kob- 
Andersen model and the Tetris model, that the so-called 
Edwards' ratio (see below for its precise definition) coin- 
cides on a wide range of densities with the Fluctuation- 
Dissipation Ratio (FDR) in homogeneous systems, i.e. in 
systems without any preferential direction. This paper 
extends those results providing a series of new evidences 
for the validity of Edwards' approach in two main direc- 
tions, (i) First of all we focus on more realistic situations 
by considering the case of granular packings subject to 
gravity: this is an important example to test the role of 



large scale inhomogeneities, such as the density profiles 
along the preferential direction, whose treatment has to 
be performed very carefully in order to avoid apparently 
non-physical results |j(J. (ii) We present new results 
concerning the independence of the FD ratio of the ob- 
servables used for its definition. 



It should be noted that previous measurements of FDR 
in the presence of gravity (and thus heterogeneities) have 
been attempted in [|o[ M. However, the negative re- 
sponse functions observed in p0| was subsequently shown 
in j|, 1 1 to be linked to memory effects , and not to 
Edwards' measure. Moreover, the measures of IBQ, I31J1 



were flawed by an incorrect definition of the correlation 
part of the Fluctuation-Dissipation ratio, due to the fact 
that one-time quantities are still evolving (see section ^ 
for a detailed discussion). The conclusion of |3l[] about 
the existence of a dynamical temperature was thus pre- 
mature. 



The link established between Edwards' approach and 
the Fluctuation-Dissipation relations could open new 
perspectives from two different points of view. First 
of all from the experimental point of view where pos- 
sibilities are open to (i) check the Edwards' measure by 
means of dynamical measurements; (ii) perform dynam- 
ical measurements (through the Fluctuation-Dissipation 
relations) of a "temperature" which should only depend 
on the density. This could allow for an at least partial 
equilibration of the disproportion existing between the 
huge number of theoretical/numerical works (and this 
paper contributes to this number) and the few experi- 
mental results. Moreover very focused experimental re- 
sults could help in discriminating between the different 
models proposed in literature. On the other hand, from 
the theoretical point of view one is left with several ques- 
tions: why Edwards' measure seems to be correct in a 
wide range of situations? Is it possible to identify some 
first principles justifications or derivations for it? Why 
the outcomes of the Edwards' approach seem to coin- 
cide with the results of Fluctuation-Dissipation measure- 
ments? 

This paper, far from being able to address all these 
questions, tries to make the link between Edwards' ap- 
proach and the Fluctuation-Dissipation measurements 
firmer in several realistic situations and propose some 
possible experimental paths for its checking. The outline 
of the paper is as follows: we first recall in section II 
the definition of the models under consideration, and 
in section III how to construct Edwards' measure. The 
case of homogeneous compaction for the Tetris model is 
described in section IV, while fluctuation-dissipation ra- 



tios during a gravity-driven compaction are measured for 
both KA and Tetris models in section M. Finally, possible 
experiments are proposed in section W% and conclusions 
are drawn. 
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II. MODELS DEFINITION 

The models we consider are lattice models, and in this 
sense are not realistic microscopic models of granular ma- 
terials. However, they are worth investigating: on the 
one hand, they have been shown to reproduce the com- 
plex phenomenology of granular media (see j?], 11, [ll], plf| 
and Q H^]); on the other hand, the validity of Edwards' 
measure for some observables has already been shown 
in an ideal case of homogeneous compaction |2^, pi[ , 
making these finite-dimensional models good candidates 
for further investigations in more realistic situations, i.e. 
with heterogeneities induced by gravity and the existence 
of a preferential direction. 



A. Tetris model 

Under the denomination of "tetris" falls a class of lat- 
tice models whose basic ingredient is the geometrical 
frustration. The models are defined on a two-dimensional 
square lattice with particles of randomly chosen shapes 
and sizes. The only constraint in the system is that par- 
ticles cannot overlap: for two nearest-neighbor particles, 
the sum of the arms oriented along the bond connecting 
the two particles has to be smaller than the bond length. 
The interactions are hence not spatially quenched but 
determined in a self-consistent way by local particle con- 
figurations. 

In the version we use (see (24|), the particles have a 
"T" -shape (three arms of length |c£, where d sets the 
bond size on the square lattice). The maximum density 
allowed is then p max = 2/3. 



B. Kob-Andersen model 

The other model we consider is the so-called Kob- 
Andersen (KA) model pq ], first studied in the context 
of Mode-Coupling theories |32| as a finite dimensional 
model exhibiting a divergence of the relaxation time at a 
finite value of the control parameter (here the density); 
this divergence is due to the presence in this model of 
the formation of "cages" around particles at high den- 
sity (the model was indeed devised to reproduce the cage 
effect existing in super-cooled liquids). 

Though very schematic, it has then been shown to re- 
produce rather well several aspects of glasses |33|, like 
the aging behaviour with violation of FDT (34], and of 
granular compaction [ |35| . 

The successful comparison of aging dynamics and pre- 
dictions of Edwards' measure was moreover shown for 
the first time for this model, in ]2J| in the idealized 
case of homogeneous compaction. On the other hand, a 
study of the violation of the FDT during the compaction 
process (under gravity) was undertaken in pjj |, and the 
existence of a dynamical temperature was advocated p5| . 



The model is defined as a lattice gas on a three dimen- 
sional lattice, with at most one particle per site. The 
dynamical rule is as follows: a particle can move to a 
neighboring empty site, only if it has strictly less than v 
neighbours in the initial and in the final position. 

Following |2Ej| , we take v = 5: this ensures that the 
system is still ergodic at low densities, while displaying a 
sharp increase in relaxation times at a density well below 
1. The dynamic rule guarantees that the equilibrium 
distribution is trivially simple since all the configurations 
of a given density are equally probable: the Hamiltonian 
is just since no static interactions exist. 

Moreover, it is also easy to consider a mixture of two 
types of particles, by considering particles of type 1 with 
a certain value v\ for the dynamical constraint, and par- 



ticles of type 2 with v% ^ v\ |35 



III. EDWARDS' MEASURE 

Edwards' approach is based on a flat sampling over 
all the blocked configurations, i.e. configurations with 
all particles unable to move. This definition therefore 
depends on the model and of the type of dynamics. For 
example, a particle is more easily blocked in presence of 
an imposed drift, e.g. gravity. 

This approach, based on the idea of describing granular 
material with a small number of parameters, leads to the 
introduction of an entropy SecIwi given by the logarithm 
of the number of blocked configurations of given vol- 
ume, energy, etc., and its corresponding density SEdw = 
SEdw/N . Associated with this entropy are the state vari- 
ables such as 'compactivity' X^ w = ^SEdw(V) and 



'temperature' T F \ 



Edw 



c)E 



~>Edw 



(E). 



The explicit construction of Edwards' measure, as well 
as of the equilibrium measure, has been described in de- 
tail in [^3[ [24] for the Tetris and the KA models. In par- 
ticular, Edwards' measure is obtained with an annealing 
procedure at fixed density. In order to select only the sub- 
set of configurations contributing to the Edwards' mea- 
sure we introduce an auxiliary temperature T aux (and 
the corresponding (3 aux — 1/T aux ) and, associated to it, 
an auxiliary energy E aux which, for each configuration, 
is equal to the number of mobile particles. A particle 
is defined as mobile if it can be moved according to the 
dynamic rules of the original model. 

In particular one measures E aux (f3 aux , p), i.e. the de- 
crease of the auxiliary energy at fixed density, perform- 
ing an annealing procedure increasing progressively (3 aux . 
From this measure one can compute the Edwards' en- 
tropy density defined by: 

SEdw(p) = s aux ((3 aU x = 00, p) = 

/>oo 

Sequil^P) / &aux(ftaux i P*)df3 aux (I) 



where e aux {(3 aux , p) is the auxiliary Edwards' energy den- 
sity and s aux ((3 aux = 0, p) = s equi i(p) is the equilibrium 
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entropy of the model. 

For the KA model, the equilibrium entropy is sim- 
ply the entropy of a lattice gas. It is worth however 
recalling that, for the Tetris model (and in general when 
the equilibrium measure is not known analytically), the 
equilibrium measure can also be obtained with an an- 
nealing procedure: also in this case one introduces an 
auxiliary temperature T' aux — 1/ f3' aux associated with an 
auxiliary energy E' aux defined as the total particle over- 
laps existing in a certain configuration. For each value 
of T' aux one allows the configurations with a probabil- 
ity given by e^~Pa^ E a^) , Starting with a large tem- 
perature T' aux one samples the allowed configurations by 
progressively decreasing T' aux . As T' aux is reduced E' aux 
decreases and only at T' aux — (no violation of con- 
straints allowed) the auxiliary energy is precisely zero. 
The exploration of the configuration space can be per- 
formed working at constant density by interchanging the 
positions of couples of particles. This procedure is used 
to compute E' aux ((3' aux , p) and e' aux (/3' aux , p) (energy per 
particle), from which one can compute the equilibrium 
entropy per particle by the expression: 



s equil(p) — s equil(Paux ~ 00 1 P) ~ 

poo 

= Sequil (flaux = 0' P) ~ / e aux(Paux J P)^Paux 

Jo 

where e' aux (f3' aux , p) is the auxiliary energy per particle. 
For the choice made for the particles one has 



s equi i(/3' aux = 0, p) = -plnp- (1 -p)ln(l-p) 
+pln4, 



(2) 



which is easily obtained by counting the number of ways 
in which one can arrange pi? particles of four different 
types on L? sites. 
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Edwards' and equilibrium entropies, computed as a 
function of density, are reported in Fig. |l| for the Tetris 
model. It is also possible to compute the so-called Ed- 
wards' ratio, defined as 



X 



Edw — 



dsEdw(p) I ds equi i{p) 



dp 



dp 



(3) 



X-Edw approaches 1 as p — > p m axi since at the maximum 
density all configurations become blocked and therefore 
equilibrium and Edwards' entropies become equivalent. 



IV. FLUCTUATION-DISSIPATION RATIO FOR 
THE CASE WITHOUT GRAVITY 

A. Tetris model: dynamics without gravity 

It is worth to recall that for the KA, the compaction 
dynamics was obtained by means of a 'piston', i.e. by cre- 
ating and destroying particles only on the topmost layer 
(of a cubic lattice of linear size L) with a chemical poten- 
tial p The validity of Edwards' measure in this case 
has been described in 
here on the Tetris mode 




FIG. 1: s Edm {p) and s equi i(p) for the Tetris model. 



We will therefore focus 
In p4| , p6| , compaction dynam- 
ics without gravity has been implemented for this model 
in order to avoid generating a preferential direction. 

The system is initialized through a random sequen- 
tial adsorption of "T" -shaped particles on an initially 
empty lattice. The grains must satisfy the geometrical 
constraints with their nearest neighbours and cannot dif- 
fuse on the lattice. This filling procedure stops when no 
other particles can be deposited on the system anymore, 
yielding a reproducible initial density of p « 0.547. 

The irreversible compaction dynamics is then realized 
alternating attempted random diffusions (in which a par- 
ticle is chosen with uniform probability and allowed to 
move toward one of its nearest-neighbours with probabil- 
ity p = j , only if all the possible geometrical constraints 
are satisfied), and attempted random depositions on the 
lattice (an empty site is chosen with uniform probabil- 
ity and a grain is then adsorbed on the lattice only if 
no violation of the geometrical constraints occurs). The 
global density increases, the system remaining homoge- 
neous during the process. 

The Monte Carlo time unit is defined as the number 
of elementary dynamical steps normalized to the num- 
ber of sites of the lattice, L 2 . In order to overcome the 
problem related to the simulation of very slow processes 
and obtain a reasonable number of different realizations 
to produce clean data, we have devised a fast algorithm 
(in the spirit of Bortz-Kalos-Lebowitz algorithm p6|), 
where the essential ingredient is the updating of a list 
of mobile particles (whose number is n rno b). In order 
to reduce the number of less significant events, such as 
failed attempts of deposition/diffusion, this algorithm is 
essentially based on a guided dynamics where only mo- 
bile particles (i.e. grains which could diffuse toward a 
neighbouring site) are considered. At each time step one 
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mobile particle is chosen with uniform probability and al- 
lowed to move if all the geometrical constraints are satis- 
fied. If the attempt has been successful, the list of mobile 
particles is then updated, performing a local control of 
grains' mobility. This procedure therefore introduces a 
temporal bias in the evolution of the system, which has 
to be taken into account by incrementing the time of an 
amount At = X/n mo b, after each guided elementary step. 
This algorithm becomes very efficient as the density of 
the system increases, since the number of mobile parti- 
cles reduces drastically. 

During the compaction, we measure the density p(t) 
of particles, the density p m ob{t) of mobile particles, the 
mobility X (t w ,t w +t) = ^J2* Ef=x ^W^+g-i^))) 
obtained by the application of a random force 
to the particles between t w and t w + t, and 
the mean square displacement B(t w ,t w + t) = 

Jn EaEti (K{t w + 1) - rt(O) 2 ) (N is the number 
of particles; a = 1, ■ • ■ ,d runs over the spatial dimen- 
sions: d = 2 for Tetris, d = 3 for KA; r£ is measured in 
units of the bond size d of the square lattice). Indeed, the 
quantities x(tw,tw +t) and B(t Wl t w + t), at equilibrium, 
are linearly related (and actually depend only on t since 
time-translation invariance holds) by 

2X(*) = is(t), (4) 

where X is the so-called Fluctuation-Dissipation ratio 
(FDR) which is unitary in equilibrium. Any deviations 
from this linear law signals a violation of the Fluctuation- 
Dissipation Theorem (FDT). Nevertheless it has been 
shown, first in mean- field models |37j , then in various nu- 
merical simulations of finite dimensional models |38|, |3^] 
how in several aging systems violations from (Q) reduce 
to the occurrence of two regimes: a quasi-equilibrium 
regime with X = 1 (and time-translation invariance) 
for "short" time separations (t <C t w ), and the aging 
regime with a constant X < 1 for large time separations. 
This second slope is typically referred to as a dynam- 
ical temperature Td yn > T^ q such that X — Xd yn — 

We have simulated lattices of linear size L = 
50, 100, 200, in order to ensure that finite-size effects 
were irrelevant. We have chosen periodic boundary con- 
ditions on the lattice, having checked that other types 
of boundary conditions (e.g. closed ones) gave the same 
results. We have investigated the irreversible compaction 
dynamics of the system up to times of 2 x 10 5 MC 
steps, realizing a large number of different runs (N runs ~ 
8000 -7- 9000) , in order to obtain clean data. The random 
perturbation is realized by varying the diffusion proba- 
bility of each particle from the initial value p = \ to the 
value p e — j + fl ■ e, where f[ — ±1 is a random variable 
associated to each grain independently for each possible 
direction (r = x,y), and e represents the perturbation 
strength. The results presented here are obtained with a 
perturbation strength e = 0.005, having checked that for 
0.002 < e < 0.01 non-linear effects are absent. 



B. Interrupted aging regime 

When the compaction process is stopped at a certain 
time t w , the system relaxes toward equilibrium: the mean 
square displacement and the integrated response function 
satisfy the time translational invariance (TTI), as Fig. ^ 
shows. This is the so-called regime of interrupted aging, 
characterized by an increase of p m ob(t) toward its equi- 
librium value (inset of Fig. ||) and a single linear relation 
for the x vs - B parametric plot. Fluctuation-Dissipation 
Theorem is then recovered, and the value of the equilib- 
rium fluctuation-dissipation ratio, Tj 9 , is obtained from 
the measures of those quantities; its value actually does 
not depend on the density of the system. 




B(t+VJ 

FIG. 2: Einstein relation in the Tetris model with interrupted 
aging, at various densities. Insets: Density of mobile particles 
vs. t, for t w = 10 4 , and B(t w + t,t w ) vs. t for two different 
t w (illustration of TTI). 

It is interesting to mention that the behavior described 
in this subsection can actually also be observed for the 
KA model, which was studied in p3| , pi| only under con- 
tinuous compaction. 



C. Aging dynamics 

If the compaction is not stopped at time t w , the den- 
sity increases, and a typical aging behaviour is observed, 
as shown in Fig. |^ where different curves of B(t + t w , t w ), 
for different values of t w , are reported: the mean square 
displacement depends explicitly not only on the observa- 
tion time t but also on t w . The system remains out of 
equilibrium on all the observed time scale, with the viola- 
tion of time translational invariance (TTI). The inset of 
Fig. H shows the collapse of the curves of B(t + t Wl t w ) for 
different values of the age of the system, obtained with 
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the following scaling function 



B(t + t w ,t w ) = c ■ 



ln( 



) 



ln(^) 



(5) 



where c, t s and r are fit parameters. Our results show a 
linear dependence of the parameters t s and r on the age 
of the system, while the coefficient c is nearly constant. 
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, tu,) for various fcu, (from bottom 
to top, 5 x 10' 3 , 10 4 , 3 x 10* and 5 x 10 4 ); inset: collapse of 
B; for the collapse, we have not taken into account the final 
portion of each curve, because of a saturation of B due to 
finite-size effects. 
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Moreover, we observe the density of mobile particles 
Pmob(t) getting smaller than the corresponding value at 
equilibrium ||24f , as another evidence of the out of equi- 
librium behaviour of the system during the compaction 
process. 

The system also features a violation of the Fluctuation- 
Dissipation Theorem. More precisely, the x vs - B para- 
metric plot, reported in Fig. shows two different linear 
behaviours: for times t smaller than t w we observe a first 
quasi-equilibrium regime where FDT holds (i.e. Xd yn = 
1), followed by a second regime in which FDT breaks 
down and a dynamical temperature Td yn arises which 
is independent of the observation time t. This quantity 
actually depends on the age of the system, t w , and there- 
fore on the density. We have investigated this behaviour 
in a large range of waiting times, corresponding to sev- 
eral values of the density, obtaining the following results: 
Xlyn = 0-646 ± 0.002, X 2 dyn = 0.767 ± 0.005, X\ yn = 
0.784 ±0.005 at = 10 4 , t 2 w = 3 x 10 4 and t 3 w = 5 x 10 4 . 

D. Results for a bi-disperse system 

In order to investigate the dependence of the dynami- 
cal temperature on the observables considered for its def- 



FIG. 4: Einstein relation in the Tetris model: plot of the 
mobility 2T^ q x(tw,t w +t) vs. the mean-square displacement 
B(t w ,t w +t), for (a) tl = 10 4 and (b) ti = 3 X 10 4 . The 
value of T| 9 is taken from Fig. ^, and the slope of the full 
straight line is one. The dashed lines are linear of the FDT 
violation whose slope gives a measure of Xdyn- 



inition, we have introduced two different types of parti- 
cles. Besides the "T" -shaped particles, already seen at 
the beginning of the section, we have considered "L"- 
shaped particles. Such grains are characterized by a 
different degree of disorder, so we expect these parti- 
cles, the "smaller" ones, to move more easily, being less 
constrained. Although we have two different types of 
particles, the system remains homogeneous during com- 
paction. 

We have therefore measured FDR for the two different 
types of particles. The results obtained show that the 
two dynamical temperatures are equal within the error- 
bars, even though the related diffusivities are different 
(see Fig. || for an example with equal fractions of T- 
shaped and L-shaped particles). This result has several 
important consequences. First of all the coincidence of 
the results for the dynamical temperature obtained with 
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different observables is a crucial step for the establish- 
ment of a thermodynamical interpretation. Another im- 
portant consequence arises from the experimental point 
of view: since the value of FDR is independent from the 
shape of the particles, it could be possible to measure 
the dynamical temperature of a granular material using 
a tracer particle different from the particles composing 
the system. This dynamical temperature should not de- 
pend on the shape of the tracer particle. The extension 
of this result to the case of compaction with gravity will 
be discussed in the next section. 




5 10 15 
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FIG. 5: Plot of the mobility 2T^ q x(tw,t w + t) vs. the mean- 
square displacement B(t w , t w + t), for different types of par- 
ticles (with equal concentrations) : (a) "T" -shaped particles 
and (b) "L"-shaped particles, t w = 5 x 10 4 . The dynamical 
temperatures associated to the two different types of particles 
(i.e. the slopes of the dashed straight lines) are equal within 
the error-bars. 



E. Comparison with Edwards' measure 



We are now in a position to compare the results for the 
Fluctuation-Dissipation Ratio Xd yn = T^ 9 /Td yn , mea- 
sured during the compaction dynamics, with the out- 
comes of Edwards' approach at the corresponding den- 
sities. In Fig. H we report the values of XecLw vs. p 
(as obtained with the Edwards' measure) and Xd yn (ob- 
tained by the FDR in dynamical measurements) at three 
different values of t w . In order to check the matching be- 
tween XEdw and Xdyn it is enough to compare the den- 
sities obtained from Fig. ^| by imposing XEdw = Xd yn 
with the corresponding dynamical densities obtained at 
the corresponding t w . From Fig. ^| one gets: p\ « 0.596 
for t 1 = 10 4 , p 2 w 0.603 for t% = 3 X 10 4 , p 3 « 0.605 
for = 5 x 10 4 . On the other hand, the evolution of 
the density of the system during the measurements of the 
FDR is reported in Fig. 0. Since the measurements are 
performed during the compaction, the density is evolv- 
ing, going from p{t w ) to p(t w + t max ). In each case, we 
obtain that indeed pi € [p(t l w ) , p(t l w + t max )], where we 
have denoted with pi the densities obtained from Fig. ^ 
for different values of [i = 1, 2, 3). More precisely, p, 
is very close to p(t l w +t max ). This is to be expected since 
the measure of the FDT violation is made for times much 
larger than t l w and, since the compaction is logarithmic, 
the system actually spends more time at densities close 
to p(t l w + t max ) than to p(t % w ). 
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FIG. 6: Static ratio XecLw as a function of density. The hor- 
izontal lines correspond to the dynamical ratios Xdyn- mea- 
sured at t w — 10 4 , 3 x 10 4 , 10 5 and determine the values 
pi w 0.596, p2 ~ 0.603, p3 w 0.605, to be compared with 
Fig. 0. 
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FIG. 7: Evolution of the density during the measurements of 
X and B, for t w = 10 4 , 3 x 10 4 , 10 5 . The evolution during 
the quasi-equilibrium part is plotted with lines, and during 
the violation of FDT with symbols. The horizontal lines cor- 
respond to the densities pi, pi, pz from fig. |^. 



V. FLUCTUATION-DISSIPATION RATIO FOR 
THE CASE WITH GRAVITY 

A. Compaction dynamics 

While the use of a compaction without gravity is useful 
to study an idealized context, real compaction due to 
shaking occurs because of gravity. 

The standard way of simulating the effect of gravity 
in a lattice model is let the particles diffuse on a tilted 
(square or cubic) lattice, with probabilities p up (resp. 
Pdown — 1 — Pup) to go up (resp. down), respecting the 
geometrical or kinetic constraints. A closed boundary is 
situated at the bottom of the simulation box, of hori- 
zontal linear size L and vertical size L z >> L (lateral 
boundary conditions can be closed or open, and vari- 
ous aspect ratios can be used, without changing quali- 
tatively the results). The control parameter is the ratio 

Pup / Pdown ^ 1 • 

This corresponds in fact to the dynamics at tempera- 
ture T = — 1/ ln(x) for the Hamiltonian 



H 



(6) 



where the Zi are the heights of the particles (of either of 
the models under consideration) above the bottom. In- 
deed, attempted diffusion moves which respect the con- 



straints are accepted with probability min(l, 



As a result, particles tend to diffuse more easily toward 
the bottom. A non-zero value of x is however needed in 
order to allow for rearrangements. 

A simple lattice gas (with the only constraint of sin- 
gle occupancy) diffusing with the above rule displays an 



equilibrium behaviour, with the known density profile 

1 



P(*) = 



I + exp(/3(z - z Q )) 



(7) 



(with [3 = 1/T and zo depends on the number of parti- 
cles) and time translation invariance, and FDT is obeyed. 
The system is therefore stationary and no evolution of the 
density occurs. 

On the other hand, systems of constrained particles 
like the Tetris or Kob-Andersen models are unable to 
reach this stationary state and are stuck at lower densi- 
ties (larger potential energies), with slow compaction and 
aging, reproducing the phenomenology of granular com- 
paction [|, [Ll| H|, j35fl . In particular, it has been shown, 
both experimentally and numerically, that, due to het- 
erogeneities, the value of the potential energy (or of the 
bulk density) is not the only relevant parameter jl0|, [ll| , 
and that, in order to explain for instance memory phe- 
nomenon, it is necessary to take into account the whole 
density profile along the vertical direction. For example, 
various density profiles can be obtained at approximately 
the same potential energy, by varying the evolution of the 
forcing T with time. 

A first order treatment consists in separating a slowly 
compacting bulk part and an interface Because the 
interface is much more dilute than the bulk, the particles 
feel much less the constraints, and it turns out that the 
density profile at the interface is exactly the same than 
the density profile of a lattice gas without constraints (see 
the inset of Fig. ||), with Hamiltonian (^J) and forcing T. 
This part of the system can therefore a priori be consid- 
ered as "in equilibrium" , i.e. its shape and dynamics are 
simply linked to the forcing. 

At this level of treatment, the system is therefore con- 
sidered to be homogeneous in the horizontal directions, 
and heterogeneities are taken into account only in the 
vertical direction. 

Compaction data, under the effect of gravity, for var- 
ious types of Tetris model can be found in g, [n], ^ 



and for the KA model in 35 . Although we have moni 



tored the usual quantities describing the compaction, we 
will therefore not repeat this analysis, but concentrate 
on the violation of FDT during compaction, showing en 
passant that the few existing data can be misleading or 
misinterpreted. 

The existence of heterogeneities along the vertical di- 
rection moreover leads to the following remarks: 

• while calculating the Edwards' measure, imposing 
only the potential energy of the system will lead 
to a unique density profile. Since the dynamical 
density profile depends on the history, it is already 
clear that specifying only the energy will not be 
sufficient to predict all dynamical observables. 

• heterogeneities exist only along the vertical direc- 
tion, so that observables along the vertical and hor- 
izontal directions should a priori be treated sepa- 
rately. 
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FIG. 8: Typical density profile obtained during the com- 
paction dynamics of the KA model. The inset shows the 
dynamical interface for x = 0.1 (circles), x — 0.2 (squares), 
x — 0.3 (diamonds), and x — 0.4 (crosses) together with equi- 
librium interfaces of a lattice gas without constraints under 
gravity at the same shaking amplitudes (lines). The interface 
profiles have been horizontally shifted for clarity. 



• dynamical measures can either be made over the 
whole system or restricted to the bulk. In the first 
case, the interface will obviously give an "equilib- 
rium" contribution that may be much larger than 
the bulk contribution. 

Our numerical results have been performed in the fol- 
lowing conditions. 

• Tetris model: Horizontal size L = 40 Number of 
particles N part = 1600 number of runs N runs — 
2000. 

• Kob- Andersen model: Horizontal size L = 20 — 40, 
Number of particles Npart — 30000—50000, number 
of runs N runs = 20 — 50. 



B. FDT and its violations in the existence of a 
preferred direction 

In general, the fluctuation-dissipation theorem relates, 
for a system at equilibrium, conjugated response and cor- 
relation functions: 



TR(t,t') 



dC(t,f) 
dt> 



Integrated between t w and t + t w in order to use the 
integrated response function, the relation becomes 

Tx{t + t w ,t w )=C(t + t w ,t + t w )-C(t + t w ,t w ) . (8) 

Let us first consider the case of horizontal degrees of 
freedom: in these directions the system is homogeneous, 



and without drift. Then if we define 
1 N 

Ch (t+t w ,t w ) = — y] (aj(t w + t)aj(t w )) 




i N 



2N ^ ^ 

a — x , y i—1 



t) (9) 



((ai(t w +t) - a 2 (t lu )) 2 ) , 



(10) 

2(C h (t + t w ,t + 



it is easily seen that Bh(t + t w ,t w ) 

t w ) Cfi (t -\- t w , tu))) ■ 

Moreover, to measure susceptibilities, a perturbation is 
applied in the following way: until t w , the system evolves 
with forcing x and Hamiltonian (||) ; at t w , a copy is made 
and evolves after t w according to the perturbed Hamil- 
tonian H e = < + H h € where H h t = + 9iVl), 

with fi, Qi = ±1 randomly for each particle and x\, y[, z\ 
are the positions of the particles in the perturbed system. 
The integrated response 



1 



JV 



Xh{t w + t,t w ) = ^^Z^ifi +*) - Si (*«;)] 

i—1 



+g l [y r l (t w + t)- y% (t w )]) (11) 

can then be measured. 

For a system at equilibrium (for example the simple 
lattice-gas with single occupancy and no kinetic con- 
straints, and Hamiltonian (||)), the FDT relation can be 
observed: 



Bh(t w +t,t w ) — 2Txh(t w +t,t w ) 



(12) 



During compaction, the violation of FDT can then be 
investigated from a parametric plot of Xh vs. B^. This 



is exactly similar to the homogeneous case of section IV 



Up to now however, the only tentative measures of 
FDR have been realized with observables coupled to the 
vertical direction plj . This is in contrast with other 
cases of systems with a preferential direction, where mea- 
sures along the only direction with no a priori hetero- 
geneities were undertaken |2j| |il], Q . 

• In |n|, the case of the KA model with a vertical 
random perturbation was considered. The vertical 
mean square displacement 

1 - 

B v (t + t w ,t w ) = —J2((z l (t w + t)-z l (t w )) 2 ) 

i=l 

was measured and confronted to the integrated re- 
sponse 

j N 

Xv(t w + t, t w ) = — ^ fi [zl(t w + t) - Zi(t w )} ^ 
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to a perturbation H% — eJ2ifi z i (/« = ^ ran ~ 
domly) . The existence of a dynamical temperature 
was inferred from the observed linear relation be- 
tween B v and \ v , with a slope different from the 
applied temperature. 

• In j30|, a perturbation in the forcing was applied, 
and confronted to the following mean-square dis- 
placement: 

B v (t + t w , t w ) = ((h(t w +t)~ h(t w )) 2 ) 

with h(t) — J2i z i(t)/N- The perturbation in the 
forcing lead to the observation of negative response 
functions (x v (t + t w ,t w ) = h r (t + t w ) - h(t + t w ), 
where h r is the mean height of the perturbed sys- 
tem, the perturbation being applied after t w ), in- 
terpreted as the signature of a "negative dynamical 
temperature" . This case was investigated in [|[ [ll[] 
where this result was shown to be linked to the 
existence of memory effects, as also confirmed in 
experiments ]h| . 

In both cases however, the existence of a downward 
drift, due to compaction, was not taken into account 
for the correct definition of the correlation part of the 
fluctuation-dissipation relation: indeed, in the first case, 
the correlation being 

1 N 

1=1 

B v (t + t w ,t w ) is not proportional to C v (t + t w ,t + t w ) — 
C v {t + t w ,t w ) as in the homogeneous case. This is even 
more easily seen in the second case, where the correlation 
conjugated to the response to a change in the driving is 
C v {t + t w ,t w ) = (h(t + t w )h{t w )) - (h(t + t w )){h{t w )). 
Indeed 

B v (t + t w ,t w ) = {h 2 (t + t w ) + h 2 (t w ))-2(h{t + t w )h(t w )) , 
and 

Cv tyj: i ~t~ ^"w) C v (t -\- t w . t w *) — 

(h 2 (t + t w )) - (h(t + t w )h{t w )) 
- (h(t + t w ))(h(t + t w )-h(t w )) (14) 

are not simply related since (h(t + t w )) ^ {h(t w )) and 
{h 2 (t + t w )) ^ (h 2 (t w )) (see also a similar discussion, on 
the case of one-time quantities changing with time, in 

It turns therefore out that the results of |3C|, |31j] are a 
priori flawed from an incorrect measure of the correlation 
part of FDR. 

We will see in the next subsections how measures of 
correlation and response functions along the horizontal 
directions lead to sensible results, whereas all measures of 
vertical correlations or response lead to the impossibility 
of defining effective temperatures. 



C. FDR in the aging (compacting) regime 

1. Vertical observables? 

Two sets of response and correlation functions can a 
priori be measured: the incoherent ones (C v ,Xv) as in 
]3l| or the coherent ones (C v , \ v ) as in |30| . 

If we write 

C v \p ~1~ t w , t -\r t w ^) C v {t -\- t W: t w ^ — 
N 



= ( j^y^^zj(t + t w )(zj(t + t w ) - Zj(t w ))j 
/ 1 N 

\ i=l 

we can observe that generically zi(t w + 1) < Zi(t w ) since 
the system is compacting, so that two opposite contribu- 
tions can be distinguished in C v (t + t w , t + t w ) — C v (t + 
t w ,t w ): the particles such that Zi(t + t w ) < h(t + t w ) 
give a positive contribution, those such that Zi(t + t w ) > 
h(t + t w ) give a negative one. At short and interme- 
diate times, the particles closer to the surface move 
more than those in the bulk and therefore the nega- 
tive contribution dominates. This leads to a negative 
C v (t + t w ,t + t w ) — C v (t + t w ,t w ). At very long times 
C v (t + t w ,t + t w ) — C v (t + t w , t w ) has to become positive 
by definition, but such times may not be reachable in a 
numerical simulation. 

This peculiar behaviour comes from the fact that the 
drift is not homogeneous in the system: some regions are 
compacting more than others. Local drifts should then 
be taken into account. However this is numerically (and 
also experimentally) too difficult to measure. 

On the other hand, the coherent correlation and re- 
sponse C v , Xv can also be measured. The difficulty 
arises from the measure of the coherent correlation func- 
tion C v , of order 1/N for N particles: relatively small 
systems have to be simulated with a large number of 
realizations. For an equilibrium lattice gas without ki- 
netic constraints, FDT is then recovered: N(C v (t + 
tw>t -\- t w ^ C v (i, -\- twi^io*)*) Xvij^ tyjjt w ^. In the 
case of the compacting system, the parametric plot of 
N(C v (t + t w , t + t w ) — C v (t + t w ,t w )) versus \ v (t -|- tyj , tyj^ 
reveals a first part of slope one, which corresponds to 
the fast, equilibrium response of the interface. At larger 
times t however, the response of the bulk, which can com- 
pactify more easily if the forcing is increased, leads to a 
decrease of Xv(t + tw,t w )i which can even become neg- 
ative as observed in Q. As t w goes to oo, the bulk 
becomes so compact that its contribution goes to zero, 
and the equilibrium FDT can be recovered thanks to the 
interface contribution. Those results are summarized in 
Fig. | 
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FIG. 9: Xv(t+tw,t w ) vs. N(C v (t+t w ,t+t w )-C v (t+t w ,t w )) 
for the KA model, with L = 20, N part = 4000, N run = 500, 
x — 0.8 (circles) and x — 0.5 (squares); t w = 2 14 and 
t — 2, • ■ ■ , 2 16 . The first equilibrium part corresponds to the 
interface dynamics, while at longer times \ v decreases be- 
cause of the bulk response. The diamonds correspond to a 
simulation with no kinetic constraint and N part = 4000 par- 
ticles: only equilibrium FDT is then observed. The straight 
line has slope 1. 

The previous investigations shows that no definition of 
an effective temperature can be inferred from dynami- 
cal measures correlated with the preferred directions in 
which heterogeneities occur. 

Note that this kind of situation also arises in the study 
of effective temperatures in driven vortex lattices with 
random pinning: while an effective temperature can be 
defined and measured for degrees of freedom perpendicu- 
lar to the drive, problems are encountered when dealing 
with longitudinal observables |Q . 

We now turn to horizontal observables. 



2. Horizontal observables 

The first result is obtained by studying the whole sys- 
tem with a horizontal perturbation applied: the rela- 
tion between mean-square displacement Bh and response 
function Xh is then clearly linear, with a slope equal to 
the temperature T of the forcing. This seemingly sur- 
prising result is easily explained by the fact that both 
functions are completely dominated by the contribution 
of the interface where the particles can diffuse quite eas- 
ily, and which actually displays an equilibrium profile (see 
Fig. ||) . It seems therefore natural to restrict the study of 
the observables to the bulk part of the sample, in which 
the density is quite homogeneous when a constant forcing 
is applied (see EL |35| and Fig. ||). The sums defining 
Bh, Xh are therefore restricted to the particles that re- 
main between t w and t w + t in the bulk (defined e.g. as 



Zmin < z i < z max , with z min and z max appropriately 
chosen) . 
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FIG. 10: 2T\h vs. B h for the Tetris model (a) and the KA 
model (t w — 2 14 )(b). The straight lines have slope 1. 

Our results, summarized in Fig. [)] and |l(| are quali- 
tatively similar for both models. The clear violation of 
FDT obtained with horizontal perturbations allows for 
the measurements of the FD ratios while nothing can be 
said using the data obtained with vertical perturbations. 



D. Results for a bi-disperse system 

In the models we have considered, it is quite easy to 
implement the presence of two types of particles (this 
has already been seen in section [lV| for the Tetris model 
without gravity). 

For the Kob- Andersen model, we can simulate "small" 
and "large" particles by taking different values for the 
kinetic constraint, e.g. v\ — 5 and — 6 or V2 = 7. 
As also shown in EjjJ, partial segregation then occurs 
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because the particle with larger v are less constrained 
and can move more easily toward the bottom. 

However, as shown in Fig. [n], there exists a bulk region 
in which the density profiles for both types of particles 
are flat. It is therefore possible to measure FDR in this 
region. 



two types of particles are equal. As already noted in sec- 



tion IV, this result is important since it means that the 
FDR can in principle be measured using tracers differ- 
ent from the particles composing the granular material, 
and that the FDR should then be independent from the 
shape of the tracer. 
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FIG. 11: Density profiles for a bidisperse system of 30000 
particles with v\ — 5 and 10000 particles with v-z = 6, for a 
forcing x = 0.2, after 2 14 time steps. The less constrained 
particles have diffused more easily towards the bottom. 
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FIG. 12: KA model: 2T\h vs. Bh, measured in the homo- 
geneous bulk, for the two types of particles (symbols: more 
constrained particles with v\ — 5; lines: less constrained par- 
ticles with V2 = 6), for various forcing (O.f , 0.2 and 0.4) and 
t w — 2 14 , t = 2, • • ■ , 2 18 . The two kinds of particles display 
different Xh and B^ at a given time, but the same violation 
of FDT. 



The results, shown in Fig. [12], are quite clear: although 
the smaller particles are more mobile than the larger 
ones, and therefore diffuse more easily, the FDR for the 



E. Edwards' measure 

Edwards' measure is defined as a flat measure over all 
blocked configurations, i.e. configurations with all parti- 
cles unable to move. For a system under gravity, a parti- 
cle at height z is 'blocked' if it cannot move downwards, 
i.e. if all its neighbouring sites at z — 1 are occupied or 
if the particle cannot move toward either of these sites 
because of the geometrical or kinetic constraints. 

Implementing this difference (with respect to the case 
without gravity) into the auxiliary model of |2^, [24] is 
straightforward. 

As in ^9|, the following procedure is used: the auxil- 
iary model has total 'energy' P aux E aux +/3E p where E aux 
is as usual the number of mobile particles, (3 aux the aux- 
iliary inverse temperature, E p the potential energy and 
j3 a Lagrange multiplier. For each value of (3, an anneal- 
ing procedure is performed on f3 aux , until configurations 
with E aux — are reached. The density profiles are 
then measured, along with the value of E p . Repeating 
the procedure then yields the curve f3Edw(E p ) directly. 
The profiles at various values of E p are shown in Fig. |f3|. 
These profiles are quite different from the dynamically 
obtained profiles at similar energies. This is not sur- 
prising since they have been obtained imposing only the 
potential energy, whereas the dynamical profiles depend 
on the history and it has been shown || |f0| [□]] that E p 
is not the only relevant parameter. 

In this case, Edwards' measure, if constructed by im- 
posing only one parameter, is not able to predict dynam- 
ical observables. 

On the other hand, since dynamically the bulk den- 
sity profiles are flat, we can generate blocked configu- 
rations with homogeneous density, at various densities. 
This yields a restricted Edwards' measure; proceeding as 
in |23|, |24| we obtain Edwards' entropy at the densities 
considered and we can therefore compute 



l(p) 



XEdw — 



dp 



dSEdw(p) 

dp 



(15) 



which is shown in Fig. [l4] in the case of the KA model. 



F. Comparison and discussion 

From the dynamics on the one hand and Edwards' 
measure on the other hand two sets of data are obtained: 

• dynamical FDR at various densities, for horizon- 
tal displacements and response functions; since the 
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FIG. 13: KA model: Density profiles for Edwards' measure, 
obtained with various values of the Lagrange parameter j3 
(symbols). A dynamically obtained profile is also shown for 
comparison: it corresponds to a constant shaking x = 0.4 and 
to a potential energy similar to the case — 0.02 (stars). The 
profiles are very different. 




FIG. 14: Circles: X E 



l(p) i dsEdui(p) 



vs. p for KA 



dp i dp 

model, imposing homogeneity; squares: dynamically obtained 
FDR (for horizontal displacements and response) vs. den- 
sity measured for the corresponding dynamical profiles. In- 
set: equilibrium and Edwards' entropy densities (imposing 
homogeneity for Edwards' entropy). 



density is evolving during the measures, an uncer- 
tainty is observed; 



statically obtained X 



Edw 



Fig. [l4| shows that the agreement between both sets of 
data is very good, even for a quite large vibration x = 0.4 
or low densities 0.73. 



The theoretical results can be summarized as follows: 
in the case of a homogeneous bulk, the ratio X& of the 
horizontal dynamical temperature to the imposed tem- 
perature — l/ln(a;) only depends on the bulk density, 
and is given by X^dw^p)- Using various vibration am- 
plitudes, we have checked that Xd at various densities 
and XEdw(p) indeed coincide. Another check of the con- 
sistency of the theoretical construction can be made by 
comparing two dynamical procedures: if a certain forc- 
ing T\ is applied until t w , and then changed to T 2 , a 
dynamical temperature Tj 2 will be obtained. While T\ 2 
depends on T2, the equality 



rpl2 



X 



Ed, 



>{p{tw)) 



should be observed, at least if T 2 is not much higher 
than Ti (in this shown in [ pT| , the bulk density 

changes suddenly by a large amount after t w ). We have 
performed such measurements and checked that this is 
indeed the case (see Fig. O). 
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FIG. 15: x vs - B following a change in the forcing from x\ to 



the slope depends on the forcing applied after 
1 ' ■ ■ " /)'. showing the equality of 



xi at t w — 

t w (t = 2, • • • , 2 1S ). Inset: T 2 \ vs 
the dynamical ratios TijTdyn- 



VI. RELATIONS WITH EXPERIMENTS 

While the link between a dynamically measured tem- 
perature and a static measure is of great theoretical im- 
portance, it is experimentally impossible to sample such 
a measure. Theoretical predictions can be checked ex- 
perimentally only through purely dynamical measures. 

Preliminary results in this direction were obtained by 
the Chicago group by measuring the volume fluctua- 
tions with respect to the steady-state volume at different 
heights of the sample. Since the curves obtained at differ- 
ent heights do not coincide the authors concluded that 
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one single observable, namely the Edwards compactiv- 
ity, cannot account for the depth dependence of the fluc- 
tuations. These conclusions are in agreement with our 
results obtained in systems with a preferential direction 
where fixing one single observable in the Edwards ap- 
proach does not allow for the prediction of the dynam- 
ical observables unless one reduces the analysis to some 
homogeneous section of the system. 

Thanks to recent theoretical progresses, new experi- 
ments have been proposed in order to check the existence 
of dynamical temperatures |23|, [ll], Q , by monitoring 
mean square displacements and mobility of tagged parti- 
cles, or tracers, in sheared super-cooled liquids or foams, 
or in sheared or tapped granular media. The existence 
of a linear relation (B vs. \) could be tested for various 
shapes, masses, etc. of the tracers, in order to check that 
this relation is indeed defining a temperature. 

Our analysis of models compacting under gravity sug- 
gests also other types of experimental possibilities. First, 
only diffusivity and mobility in the direction perpendic- 
ular to the gravity should be measured. Moreover, the 
existence of strong heterogeneities implies that a tracer 
close to the interface should allow to measure a tempera- 
ture T\ which depends directly on the driving amplitude 
(in the models, T = — l/ln(x)), and is stationary, i.e. 
does not depend on the bulk density. On the other hand, 
a tracer well immersed in the bulk should yield another 
value T\ of the FDR (which depends on t w ), and the ra- 
tio Tj/Tj should depend only on the density of the bulk. 
For example, if experiments are performed changing the 
driving intensity at t w from x\ to and the two associ- 
ated "interface" and "bulk" temperatures are measured, 
the ratio T\jT l d should be independent of x 2 provided 
the bulk density does not change significantly. In this 
way, the comparison of only dynamical measures would 
be a strong experimental test for the whole theoretical 
construction, without any need to sample the underlying 
static measure. 



VII. SUMMARY AND CONCLUSIONS 

In this paper, we have studied two paradigmatic mod- 
els for the compaction of granular media. These mod- 



els consider particles diffusing on a lattice, with either 
geometrical or dynamical constraints. Idealized com- 
paction without gravity has been implemented for the 
Tetris model, and compaction with a preferential direc- 
tion imposed by gravity has been studied for both mod- 
els. The possibility to define dynamically a tempera- 
ture in the framework of Fluctuation-Dissipation rela- 
tions and to link it to the statically constructed Edwards' 
measure has been investigated. 

In the first, ideal, case of the homogeneous compaction, 
the obtained FD ratio have been clearly shown to be in 
agreement with the prediction of Edwards' measure at 
various densities. 

The situation is more complicated if a preferential di- 
rection is present: then the whole density profile has a 
priori to be taken into account. Moreover, the vertical 
drift due to compaction leads to contradictory (and some- 
times meaningless) results when observables coupled to 
the preferential direction are considered for the evalua- 
tion of a FD ratio 30 , 3lJ . Since the energy of the system 
is not the only parameter, and since the density profiles 
depend on the history, Edwards' measure is not a priori 
able to predict the dynamical configurations. If how- 
ever the homogeneity of the bulk is imposed, FD ratio 
obtained dynamically for horizontal displacements and 
mobility and from Edwards' measure coincide. 

It is striking to note that Edwards' measure, which 
a priori could be valid only for very weak forcing and 
almost stationary systems seems however to yield good 
predictions even for non-stationary systems that are still 
compacting. 

Finally, we have proposed experimental tests of the 
whole theoretical construction, through the comparison 
of various types of dynamical measurements, since the 
construction of Edwards' measure, numerically straight- 
forward, is obviously impossible in experiments. Ac- 
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